home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / chbev.z / chbev
Encoding:
Text File  |  2002-10-03  |  4.7 KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCHHHHBBBBEEEEVVVV((((3333SSSS))))                                                            CCCCHHHHBBBBEEEEVVVV((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CHBEV - compute all the eigenvalues and, optionally, eigenvectors of a
  10.      complex Hermitian band matrix A
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE CHBEV( JOBZ, UPLO, N, KD, AB, LDAB, W, Z, LDZ, WORK, RWORK,
  14.                        INFO )
  15.  
  16.          CHARACTER     JOBZ, UPLO
  17.  
  18.          INTEGER       INFO, KD, LDAB, LDZ, N
  19.  
  20.          REAL          RWORK( * ), W( * )
  21.  
  22.          COMPLEX       AB( LDAB, * ), WORK( * ), Z( LDZ, * )
  23.  
  24. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  25.      These routines are part of the SCSL Scientific Library and can be loaded
  26.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  27.      directs the linker to use the multi-processor version of the library.
  28.  
  29.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  30.      4 bytes (32 bits). Another version of SCSL is available in which integers
  31.      are 8 bytes (64 bits).  This version allows the user access to larger
  32.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  33.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  34.      only one of the two versions; 4-byte integer and 8-byte integer library
  35.      calls cannot be mixed.
  36.  
  37. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  38.      CHBEV computes all the eigenvalues and, optionally, eigenvectors of a
  39.      complex Hermitian band matrix A.
  40.  
  41. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  42.      JOBZ    (input) CHARACTER*1
  43.              = 'N':  Compute eigenvalues only;
  44.              = 'V':  Compute eigenvalues and eigenvectors.
  45.  
  46.      UPLO    (input) CHARACTER*1
  47.              = 'U':  Upper triangle of A is stored;
  48.              = 'L':  Lower triangle of A is stored.
  49.  
  50.      N       (input) INTEGER
  51.              The order of the matrix A.  N >= 0.
  52.  
  53.      KD      (input) INTEGER
  54.              The number of superdiagonals of the matrix A if UPLO = 'U', or
  55.              the number of subdiagonals if UPLO = 'L'.  KD >= 0.
  56.  
  57.      AB      (input/output) COMPLEX array, dimension (LDAB, N)
  58.              On entry, the upper or lower triangle of the Hermitian band
  59.              matrix A, stored in the first KD+1 rows of the array.  The j-th
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCHHHHBBBBEEEEVVVV((((3333SSSS))))                                                            CCCCHHHHBBBBEEEEVVVV((((3333SSSS))))
  71.  
  72.  
  73.  
  74.              column of A is stored in the j-th column of the array AB as
  75.              follows:  if UPLO = 'U', AB(kd+1+i-j,j) = A(i,j) for max(1,j-
  76.              kd)<=i<=j; if UPLO = 'L', AB(1+i-j,j)    = A(i,j) for
  77.              j<=i<=min(n,j+kd).
  78.  
  79.              On exit, AB is overwritten by values generated during the
  80.              reduction to tridiagonal form.  If UPLO = 'U', the first
  81.              superdiagonal and the diagonal of the tridiagonal matrix T are
  82.              returned in rows KD and KD+1 of AB, and if UPLO = 'L', the
  83.              diagonal and first subdiagonal of T are returned in the first two
  84.              rows of AB.
  85.  
  86.      LDAB    (input) INTEGER
  87.              The leading dimension of the array AB.  LDAB >= KD + 1.
  88.  
  89.      W       (output) REAL array, dimension (N)
  90.              If INFO = 0, the eigenvalues in ascending order.
  91.  
  92.      Z       (output) COMPLEX array, dimension (LDZ, N)
  93.              If JOBZ = 'V', then if INFO = 0, Z contains the orthonormal
  94.              eigenvectors of the matrix A, with the i-th column of Z holding
  95.              the eigenvector associated with W(i).  If JOBZ = 'N', then Z is
  96.              not referenced.
  97.  
  98.      LDZ     (input) INTEGER
  99.              The leading dimension of the array Z.  LDZ >= 1, and if JOBZ =
  100.              'V', LDZ >= max(1,N).
  101.  
  102.      WORK    (workspace) COMPLEX array, dimension (N)
  103.  
  104.      RWORK   (workspace) REAL array, dimension (max(1,3*N-2))
  105.  
  106.      INFO    (output) INTEGER
  107.              = 0:  successful exit.
  108.              < 0:  if INFO = -i, the i-th argument had an illegal value.
  109.              > 0:  if INFO = i, the algorithm failed to converge; i off-
  110.              diagonal elements of an intermediate tridiagonal form did not
  111.              converge to zero.
  112.  
  113. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  114.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  115.  
  116.      This man page is available only online.
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.